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Abstract 

We investigate the cogrowth and distribution of geodesies in R. Thompson's 
group F. 

1 Introduction 

In this article we study the cogrowth and distribution of geodesies in Richard Thomp- 
son's group F, in an attempt to decide experimentally whether or not F is amenable. 

The cogrowth of a finitely generated group G is denned as follows. Suppose 
S = {ai, . . . ,ak} generates and consider the Cayley graph Q of (G, S). Let r n 
be the number of paths in this graph of length n starting and ending at the identity 
element — let us call such paths returns. Since we can concatenate any two such 
paths to get another we have 

r n r k < r n+k (1) 
and then by Fekete's lemma (see, for example, [26] ) 

p = lim sup r^" (2) 

n— >oo 

*The first author acknowledges support from ARC project DP110101104. The second author 
thanks NSERC of Canada for financial support. 

1 Formally, we consider G as the epimorphic image from the free monoid generated by S U S , 
rather than S U S" -1 as being a subset of G 



1 



exists. Since we consider generators and their inverses to label distinct edges in Q, 
then p < 2k. 

The connection between this growth rate and amenability was established by 
Grigorchuk and independently by Cohen: 

Theorem 1 ( [121 E])- Let G,S and p be as above. G is amenable if and only if 
p = 2k. 

Let p n be the number of returns of length n on Q which do not contain immedi- 
ate reversals. Again concatenation shows that p n is supermultiplicative so Fekete's 
lemma gives 

a = lim sup pl( n (3) 

exists. This constant is called the cogrowth for (G, S). In this case since there are 
2k(2k — l) n ~ 1 freely reduced words of length n in the 2k generators and their inverses, 
we have a < 2k — 1. 

The previous theorem can then be restated as: 

Theorem 2 ( [T2l [7]). Let G,S and a be as above. G is amenable if and only if 
a = 2k-\. 

Note that lim sups are required since, for example, if G has a presentation where 
all relators have even length, the number of returns of odd length (with or without 
immediate reversals) is 0. 

In this article we compute bounds on the cogrowth rates of a number of 2- 
generator groups: Thompson's group F, the free and free abelian groups on 2 genera- 
tors, Baumslag-Solitar groups, and various wreath products. Each of these examples, 
apart from F, is known to be either amenable or non-amenable. We compare the 
data obtained for F against these examples, to see whether F behaves more like an 
amenable or a no n- amenable group. 

The question of the amenability of Thompson's group F has captivated many 
researchers for some time, initially since F has exponential growth but no nonabelian 
free subgroups, making it a prime candidate for a counterexample to von Neumann's 
conjecture that a group is non-amenable if and only if it contains a nonabelian free 
subgroup. In 1980 Ol'Shanskii constructed a finitely generated non-amenable group 
with no nonabelian free subgroups [15], and in '82 Adyan gave further examples pQ. 
In 2002 Ol'shanksii and Sapir constructed finitely presented examples [16]. In spite 
of these results the amenability or non-amenability of F remains an intensely studied 
problem. In 2009 two articles appeared on the ArXiv, one asserting F is amenable 
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[2H] and one claiming F is non-amenable [5]. Both papers contained errors, which 

despite some effort do not appear to be fixable [251 Ell- 
in the second half of the article we extend our techniques to study the distribution 

of geodesic words in Thompson's group. 

This work is inspired by a previous paper by Burillo, Cleary and Weist [5], who 

also applied computational techniques to consider the amenability of F. We refer 

the reader to this paper for more background on Thompson's group and the problem 

of deciding its amenability computationally. 

The present article is organised as follows. In Section 2 we compute rigorous lower 

bounds on the cogrowth by computing the dominant eigenvalue of the adjacency 

matrix of truncated Cayley graphs. We then extrapolate these bounds to estimate 

the cogrowth and compare and contrast those extrapolations for F and other groups. 

In Section 3 we use a weighted random sampling of random words in the generators 

to estimate the exponential growth rate of trivial words in several different groups. 

As a bi-product we estimate the distribution of geodesic lengths as a function of 

word length. 

2 Bounding returns and cogrowth 
2.1 Bounding the number of returns 

Consider the Cayley graph Q of some group G with finite generating set — for the 
discussion at hand, let us assume that G is generated by two nontrivial elements a, b. 

As noted above, an upper bound for the rate of growth of returns p is 4. We can 
compute lower bounds for the number of returns, and subsequently the cogrowth, as 
follows. 

Consider the following sequence of finite connected subgraphs, Qn of N vertices 
that contain the identity. 

Set Qi = the identity vertex. Record the list of edges incident to G\. Define 
Qi-, Qs, • • • by appending edges from this list, one at a time. Once the list is exhausted 
(so Qm = B(l), where B(r) is the ball of radius r in Q), repeat the process. It follows 
that for each iVeN there is an R e N so that B(R) C Q N C B(R + 1). 

We can then define r^, n to be the number of returns of length n in Qn- Since 
Qn C Qn+i, the sequence {rN,n} is supermultiplicative, so pn = lim sup^^ r% n n 
exists by Fekete's lemma. Further we must have r n > r^.n and so p > pn- Hence we 
can bound p by computing p^. 

Using the Perron-Frobenius theorem (in one of its many guises — Proposition 
V.7 from [TJJ for example) the growth rate Pn °f such paths on Qn is given by the 
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dominant eigenvalue of the corresponding adjacency matrix, provided it is irreducible. 
Since Qn is connected, the corresponding adjacency matrix will be irreducible. 

In some cases we can also demonstrate that the adjacency matrix is aperiodic, 
which implies that the dominant eigenvalue is simple and dominates all other eigen- 
values. This also implies that the corresponding generating function YlPN,nZ n has a 
simple pole at the reciprocal of that eigenvalue. Unfortunately many of the matrices 
we study are not aperiodic, but they do have period 2. 

Perhaps the easiest way to prove that the matrix is aperiodic is to show the exis- 
tence of two circuits of relatively prime length (see chapter V.5 in [11] for example). 
In order to show that a matrix has period 2 it suffices (providing the matrix is finite) 
to show that if there is a path of length k between any two nodes, then there is a 
path of length k + 2£ between those two nodes for any £. 

It follows that the adjacency matrix for Qn is aperiodic whenever the group G 
has a presentation with an odd length relator (since aa~ l and the odd length relator 
form circuits of relatively prime lengths) and if all relators have even length, the 
matrix has period 2 (since a path of length k can be made into a path of length 
k + 2£ by inserting (aa -1 )^). 

In particular we have that subgraphs of Baumslag-Solitar groups BS(p, q) with 
p + q odd are aperiodic, while subgraphs of BS(p, q) with p + q even, Thompson's 
group F, ZlZ and F 2 (the free group on 2 generators), all with the usual generating 
sets, have period 2. 

Since the above groups grow exponentially, and B(R) C Q N C B(R + 1), then 
the radius of Qn is O(logiV). 

We used this method to compute pn for a selection of groups. However, we found 
significantly better bounds by considering only freely reduced words, i. e. paths that 
do not contain immediate reversals, essentially since there is less to count. 

2.2 Bounding the cogrowth 

Let p n be the number of returns of length n on Q which do not contain immediate 
reversals. We similarly define p^^ n to be similar paths on the subgraph Q N . Again 
we define the exponential growth of these quantities by 

a = limsupp^ n a N = limsupp^ n n 

n— >oo n— >oo 

and a > a^. 

In this case, we cannot now simply concatenate two freely reduced paths to obtain 
another freely reduced path since it may create an immediate reversal. Thus we do 
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not have similar supermultiplicative relations. We can, however, relate r n to p n and 
p to a using the following result of Kouksov [13J which we have specialised to the 
case of 2 generator groups. 

Lemma 3. (from [13]) Let R(z) = ^2r n z n and C(z) = YlPnZ n be the generating 
functions of returns and freely reduced returns respectively. Then, 

1 - z 2 ( z \ 
C(z) = -R I ) and equivalently 

R ( z ) = 1 C 



16z 2 \ 6z 

A careful generating function argument gives the second equation (and the first 
is simply its inverse). Consider any freely reduced returning path; it can be mapped 
to an infinite set of returning paths by replacing each edge s by any returning path 
in the free group on 2 generators that does start with s _1 . At the level of generating 
functions, this is exactly the substitution z !->■ l y,1 6 ~ 12z2 - 

A very general result for generating functions then links the dominant singularity 
of R{z) to the value of p: 

Theorem 4. ([11], page 240.) If f(z) is analytic at and p is the modulus of a 
singularity nearest to the origin, then the coefficient f n = [z n ]f(z) satisfies: 

limsup \f n \~ 1/n = P 

Combining these two results (and using the positivity of r n ,p n ) we obtain 
Corollary 5. The constants p and a are related by 

a 2 + 3 



P 



a 

/3 2 +3 



Further if [3 is a lower bound for a, then [3' = is an lower bound for p. 

We are unable to prove a similar exact relationship between p^ and a^, but we 
do have the following bound: 

Lemma 6. For a fixed value of N we have pn < 
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Proof. Consider the generating functions of returns and freely reduced returns on 

Qn- 

Rn{z) = ^2r N:n z n C N {z) = ^2p NnZ n 

n>0 n>0 

It suffices to show that 

I — z 2 ( z \ 
C N {z) > 1 + 3z2 Rn y 1 + 3z 2 J or equivalent^ 

D -1 + 2^1-12^2 ^i-Vl-12^ 2 
Rn(z) < — 2 C N 



1 — 16z 2 \ Qz 

within the respective radii of convergence. Let u be any freely reduced returning path 
of length k in Qn — this path contributes z k to the generating function Cn(z). The 
substitution z h> 1 ~ v ' 1 q~ 12z2 maps lu to an infinite set of non-reduced returning paths 
by replacing each edge with freely reduced words from F 2 . Some of the resulting 
words will lie entirely within Qn and so be enumerated by the generating function 
Rn(z). However an infinite number of these words will not be contained in Qn- These 

words are enumerated by ~ 1+ ^ 6 ~2 2z2 C N (^p^) but not by R N {z). Thus the 
inequality follows. □ 

To compute «jv we relate it to the dominant eigenvalue of an adjacency matrix. 
Unfortunately there is no simple way to reuse the adjacency matrix of Qn, m order 
to enumerate paths without immediate reversal. Instead we construct a new graph 
Hn which encodes freely reduced paths in Q as follows: Hn has iV vertices labeled 
by pairs (1, — ) or (g, s) where g G G and s G S. The vertex (1, — ) corresponds to 
being at the identity vertex of Q, and (g, s) to being at the group element g G Q after 
having just read a letter s. The edges oiV, N are 

E{U N ) = {((g, s), (h, t)) G (V(U N )) 2 \h = gt and st ^ 1} (4) 

So a path (1, — ), si), (g 2 , s 2 ), ■ ■ ■ {gk, Sfc) corresponds to a path in Q starting at 1 
with g 1 = si, g 2 = s±s 2 , . . . g^ = sis 2 . . . Sk a freely reduced word. 

We construct Hn using a breadth- first search similar to the construction of Qn, 
starting with "Hi = (1, — ) and appending vertices one at a time so that {g G 
G | (g,s) G V,n} lies between two balls of a given radius. It follows that Hn is 
necessarily connected, and the corresponding adjacency matrices are irreducible. 

We then compute the growth rate of paths (and so freely reduced returns) oyl'Hn 
by computing the dominant eigenvalue of the corresponding adjacency matrix. 
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2.3 Exact lower bounds 



Since each node of Qn has outdegree at most 4 and those of the Hn (excluding (1, — ) 
of outdegree 4) have outdegree at most 3 (vertices on the boundary may have smaller 
degree) the corresponding adjacency matrices are sparse. We found that the power 
method and Rayleigh quotients (see [2T] for example) converged very quickly to the 
dominant eigenvalue and so the growth rate. 

We constructed Qn and Hn for many different values of iV ranging between 10 2 
and 10 7 . Our calculations on Thompson's group F as well as the Baumslag-Solitar 
groups BS(2, 2), BS{2, 3) and BS(3, 5), yielded the following result. 

Theorem 7. The following are exact lower bounds on the cogrowth, a, of the indi- 
cated groups. 

BS(2, 2) > 2.5904 BS(2, 3) > 2.42579 BS(3, 5) > 2.06357 
Thompson's group > 2.17329 

This implies that the growth rate of all trivial words, p, in these groups are bounded 
as indicated. 

BS(2, 2) > 3.78522 BS(2, 3) > 3.66250 BS(3, 5) > 3.51736 
Thompson's group > 3.55368 

Note that all of these bounds were computed using information from Hn and 
Corollary [5j We observed that the bounds obtained from Qn were worse — typically 
differing in the second or third significant digit. We also note that the above result 
for BS(2, 3) is improves on a result in [Q] (the preprint was withdrawn by the authors 
since it contained an error). 

These computations were done on a desktop computer using about 4Gb of mem- 
ory. It should be noted that while our techniques require both exponential time and 
memory, it was memory that was the constraining factor. We did implement some 
simple space-saving methods. Perhaps the most effective of these was to store ele- 
ments as geodesic words in the generators rather than as their more standard normal 
forms (eg tree-pairs for Thompson's group or words in the normal form implied by 
Briton's lemma for the Baumslag-Solitar groups). These geodesic words could then 
be stored as bit-strings rather than ASCII strings. We believe that by running these 
computations on a computer with more memory we could improve the bounds, but 
the returns are certainly diminishing. 
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2.4 Extrapolation and comparison 

The results of the previous section can be extended by considering the sequence 
of lower bounds and using simple numerical methods to extrapolate them to 
iV — > oo. This required only minimal changes to our computations; after computing 
the adjacency matrix of for the maximal value of N, we computed the dominant 
eigenvalue of submatrices. The corresponding estimate of the eigenvector was then 
used as an initial vector for estimating the eigenvalue of the next submatrix. This 
meant that we could compute a sequence of lower bounds in not much more time 
than it took to compute our best bounds. 

In Figure [I] we have plotted against Yiogiv for three non-amenable Baumslag- 
Solitar groups and F. We found that this gave approximately linear plots and so 
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Figure 1: A plot of cogrowth lower bounds against 1 /\ogN. We see that the 
groups that known to be non-amenable are converging to numbers strictly below 
3. The Thompson's group sequence has a clear upward inflection (as iV — > oo or 
yiogiv — > 0) and so it is difficult to estimate whether the limit is 3 or less than 3. 

this suggests that 

« AT ~ ftoo - A /log7V. 
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Since is a monotonically increasing sequence and is bounded above by 3, we have 
that a at — > exists. Unfortunately we cannot prove that = a, but certainly 

ctoo < a. 

Note that the curves terminate at iV = 10 7 but start at different iV-values. This 
is because the graphs Qjq do not contain freely- reduced loops for small values of N. 
The smallest value of N for which contains a freely reduced loop depends on the 
length of the relations of the group and on the details of the breadth-first search used 
to construct the graph. 

One can observe that the Baumslag-Solitar groups all seem to behave similarly 
and that the sequences of bounds are clearly converging to constants strictly less 
than 3. This is completely consistent with the non-amenability of these groups. 
Thompson's group behaves quite differently — in particular we see that the curve 
has some upward inflection (as x — > 0) and it makes it very unclear as to whether or 
not converges to 3 or below 3. 

For the sake of comparison we decided to repeat the above analysis for a set 
of amenable groups and so we computed similar sequences of lower bounds for 
BS(1, 2), BS(1, 3), Z 2 and Z I Z. These results are plotted in Figure § 

Note that no sequence gives a perfectly straight line and so to estimate we 
fitted the data to the form 

a n = a (Xl + A / (log AO 4 . 

We varied the number of data points by removing small- N points and we also varied 
the value of 5. For any fixed number of points we varied 5 to find a value that 
minimsed the R 2 statistic. This gives an "optimal" value of and A. 

For some groups, we found that these optimal values were quite sensitive to 
changes in S, while other groups were quite robust. To include some measure of this 
systematic error we moved 5 through a range of values so that the R 2 statistic was 
allowed to move to 5% below its optimal value. These results are summarised in 
Tables [Hand [2 

The results for all all the groups except Thompson's group are as one might 
expect — the amenable groups all give estimates of close to 3, and the non- 
amenable groups give < 3. Hence it would appear as though this technique is a 
reasonable test to differentiate amenable and non-amenable groups. Unfortunately 
it is not sufficiently sensitive to determine the amenability of Thompson's group. In 
particular we find that the results are too sensitive to variations in 5 and to removal 
of low-iV data points. 

Because of this, we turn to numerical methods based on random sampling and 
approximate enumeration. 



9 



Group 


Number of points 


Optimal R 2 Value 


5 range 


estimate 


BS(2,2) 


4501 

oc;nn 
ZOUU 


0.998 

n oofi 
U.yyo 


1.74 ±0.05 

1 Q K -L- n 1 Q 
l.oO ± U. lo 


2.682 ±0.007 
z.o i a ± u.uuy 


BS(2,3) 


4501 
2500 


0.999 
0.999 


1.36 ±0.04 
1.57 ±0.07 


2.597 ±0.012 
2.562 ± 0.009 


BS(3,5) 


4101 
2000 


0.998 
0.998 


1.33 ±0.05 
1.65 ±0.19 


2.29 ±0.01 
2.24 ±0.03 


Thompson 


3947 
2000 
1700 


0.998 
0.998 
0.998 


0.83 ±0.07 
0.93 ±0.16 
0.65 ±0.21 


2.79 ±0.08 
2.69 ±0.12 
2.95 ±0.38 



Table 1: Results of fitting eigenvalue data for non-amenable groups. The Baumslag- 
Solitar groups all give good results, but Thompson's group does not. There is some 
upward drift in the estimate of as one cuts out small N data, but at the same 
time the error in the estimates blows up. 



Group 


Number of points 


R 2 Value 


S range 


ttoo estimate 


BS(1,2) 


4501 
2500 


0.99975 
0.99981 


1.7316 ±0.0225 
1.9472 ± 0.0552 


3.0158 ±0.0031642 
2.9975 ± 0.0031542 


BS(1,3) 


4501 
2500 


0.99894 
0.99855 


1.354 ±0.046 
1.54 ±0.151 


3.0722 ± 0.016473 
3.0261 ± 0.026664 


Z 2 


4501 
2500 


0.99613 
0.99932 


5.1624 ±0.1134 
10.996 ±0.154 


3.002 ± 0.000364 
3 ± 1.7248 x 10~ 6 


ZlZ 


3947 
2000 


0.99925 
0.99915 


0.8592 ± 0.0436 
1.0237 ±0.1251 


3.1807 ±0.050903 
3.0476 ± 0.082052 


Thompson 


3947 
2000 
1700 


0.99796 
0.99869 
0.99848 


0.83 ± 0.072 
0.9344 ±0.1548 
0.6464 ±0.2016 


2.7866 ± 0.083778 
2.6917 ±0.12318 
2.9532 ± 0.38051 



Table 2: Results of fitting eigenvalue data for amenable groups. All the amenable 
groups give good results quite close to 3, though ZlZ is not as good as the others. Also 
note that since balls in Z 2 grow quadratically with radius rather than exponentially, 
better results can be obtained by fitting against 1 /n s rather than l /{\ogN) 5 . 
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Figure 2: A plot of cogrowth lower bounds against yiogjv. We see that the 
groups that known to be amenable are clearly converging to 3. Again we see that 
Thompson's group behaves quite differently. 

2.5 An aside — cogrowth series 

As a biproduct of our computations we obtained the first few terms of the cogrowth 
series for all of these groups. It is well know that the number of trivial words in Z, 2 
is given by ( 2 ^) (see A002894 [24] ) ; the corresponding generating function is not 
algebraic and is expressable as a complete elliptic integral of the first kind. The 
number of trivial words in F 2 is just the number of returning paths in a quadtree 
and its generating function is 3(1 + 2\J\ — \2z 2 )~ l (see A035610 [24J). 

Unfortunately we have been unable to find (using tools such as GFUN [22J) 
any useful explicit or implicit expressions for the cogrowth series (or the generating 
functions) for any of the other groups we have examined. For completeness we include 
our data in Table [3] 
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n 


Thompson 


BS(1,2) 


BS(1,3) 


BS(2,2) 


BS(2,3) 


BS(3,5) 







1 


1 


1 


1 


1 


1 


1 


1 
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5 





10 

















6 








12 


12 











7 





20 








14 








8 





64 


40 


40 








16 


9 





96 








28 








10 


20 


338 


264 


224 


60 


20 


72 


11 





736 








84 








12 


64 


2052 


1604 


1236 


240 


64 


272 


13 





5208 








564 








14 


336 


13336 


9748 


7252 


1090 


280 


1504 


15 





36330 








2760 








16 


1160 


92636 


61720 


41192 


6492 


1048 


8576 


17 





248816 








13496 








18 


5896 


665196 


412072 


247272 


33728 


4660 


46080 


19 





1771756 








75768 








20 


24652 


4776094 


2750960 


1491136 


174760 


17964 


257160 


21 





12848924 








411234 








22 


117628 


34765448 


18725784 


9119452 


958364 


77508 


1475592 



Table 3: The first few terms of the cogrowth series C(z) for various groups, i.e. the 
number of freely reduced words equivalent to the identity. The first few terms of the 
returns series R(z) can be obtained from the above using Lemma [3j 
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3 Distribution of geodesic lengths 



In this section we broaden our study from the growth rate of trivial words to the 
distribution of geodesic lengths of all words by sampling random words. In previous 
work of Burillo et al [5], random words in Thompson's group F were sampled using 
simple sampling; words were grown by appending generators one-by-one uniformly 
at random. Those authors observed only very trivial words and so then sampled 
uniformly at random from a subset of those words, namely the set of words with 
balanced numbers of each generator and their inverses. Again, very few trivial words 
were observed. Indeed if Thompson's group is non-amenable, the probability of 
observing a trivial word using simple sampling will decay exponentially quickly. 

We will proceed along a similar line but using a more powerful random sampling 
method based on flat-histogram ideas used in the FlatPERM algorithm [TBI !19.J. 
Each sample word is grown in a similar manner to simple sampling — append one 
generator at a time chosen uniformly at random. The weight of a word of n symbols 
is simply 1, so that the total weight of all possible words at any given length is just 
4 n . As the word grows we keep track of its geodesic length. We now deviate from 
simple sampling by "pruning" and "enriching" the words. 

Consider a word of length n, geodesic length £ and weight W. If we have "too 
many" samples of such words, then with probability 1 /2 prune the current sample 
or otherwise continue to grow the current sample but with weight 2W. Similarly if 
we have "too few" samples of the current length and geodesic length, then enrich 
by making 2 copies of the current word and then growing a sample from both each 
with weight Of course, one is free to play around with the precise meaning of 
"too few" or "too many" . We refer the reader to [TBJ [19] for more details on the 
implementation of this algorithm. The mean weight (multiplied by 4") of all samples 
of length n and geodesic length £, c n ^ is then an estimate of the number of such 
words. 

In order to run the above algorithm we need to be able to compute the geodesic 
length of the element generated by a given random word. Computing geodesic lengths 
from a normal form is, in general, a very difficult problem and remains stubbornly 
unsolved for many interesting groups, such as BS(2, 3). Because of this we restrict 
our studies to Thompson's group and a number of different wreath products. 

• Thompson's group — a method for computing the geodesic length of an element 
from its tree-pair representation was first given by Fordham [4j, though we 
found it easier to implement the method of Belk and Brown [3] . 

• Wreath products — we use the results of [6J to find the geodesic lengths in 
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Z I Z, Z i (Z ? Z) and Z I F 2 . 



We note that polynomial time algorithms to compute geodesies in Baumslag-Solitar 
groups has recently been given in the cases BS(l,n) [TU] and BS(n, kn) [8], but we 
have not implemented these. 

3.1 Distributions 

We used the random sampling algorithm described above to estimate the distribution 
of geodesic lengths in Thompson's group F, as well as Z I Z, Z ? F 2 and Z I (Z I Z). 
Each run took approximately 1 day on a modest desktop computer. To visualise the 
results, we started by normalising the data by dividing by the total number of words 
(i.e. 4 ra or 6 n ). The resulting peak-heights still decay with length, and we found that 
multiplying by \fn compensated for this. The normalised distributions are plotted 
in Figures [3j [4] and [5j 

In each case we see similar behaviour. At short word lengths (i.e. small n) the 
distribution of geodesic lengths is quite wide, but settles to what appears to be a 
bell-shaped distribution at moderate lengths. This suggests that the geodesic length 
has an approximately Gaussian distribution about the mean length and that the 
tails of the distribution are exponentially suppressed. This also explains why the 
normalising factor of \Jn works well. 

If this is indeed the case, then we expect that trivial words, having geodesic length 
zero, will be exponentially fewer than 4™ — implying that Thompson's group is non- 
amenable. Unfortunately things cannot be so simple, because the same reasoning 
would imply that Z \ Z is non- amenable. 

One obvious difference between the graphs is the movement of the peak of the 
distribution, that is the rate of growth of the mean geodesic length. It is clear that 
the mean geodesic length of Z I F 2 grows linearly, and so the group has a non-trivial 
rate of escape — exactly as one would expect of a non-amenable group. Similarly we 
see that the mean geodesic lengths of the other wreath products grow sublinearly, 
so their rates of escape are zero. When we examine the movement of the peak of 
Thompson's group's distribution, things are less clear; the mean geodesic length 
appears to be very nearly linear. 

Estimating the mean geodesic length for Thompson's group was substantially 
easier. We constructed 2 12 random words of length 2 16 . As each word was con- 
structed generator-by-generator, the geodesic length was computed and added to 
our statistics. So while there is correlation between the geodesic lengths at different 
word lengths within a given sample, there is no correlation between samples. This 
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Figure 3: A plot of the normalised distribution of the number of words c n ^ of length 
n and geodesic length I in Thompson's group F. Notice that the peak position is 
quite stable, indicating that the mean geodesic length grows roughly linearly with 
word length. 

took approximately 3 days on a modest desktop computer. Our data is plotted in 
Figure 6. 

We assume that the mean geodesic length grows as n v . Linear regression on a 
log-log plot estimates u ~ 0.98. Further, if we fit a moving "window", we find that 
the local estimates of v increase as the positioning of the window increases. This 
strongly suggests that the mean geodesic length grows linearly. 

To test linearity further, we generated a small number words of length 2 20 = 
1048576. It took approximately 1 hour to generate each word and compute the 
corresponding geodesic length, so this was too slow to generate meaningful statistics. 
In each case we observed that the ratio t jn appeared to converge to approximately 
0.28. Of course, this does not preclude more exotic sublinear behaviour such as 
n u (logn) 6 . Such logarithmic corrections are extremely difficult to detect or rule out. 

We now estimate the rate of escape by assuming linear growth with a polynomial 
sub dominant correction term 

(£) n = An + bn 5 . (5) 
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t/n 



Figure 4: A plot of the normalised distribution of the number of words c n ^ of length 
n and geodesic length I in Z \ 7L. Observe that the peak position is clearly moving 
towards the left of the plot suggesting that the mean geodesic length grows sublin- 
early. 




Figure 5: Plot of the normalised distribution of the number of words c n / of length n 
and geodesic length i in Z \ F 2 (left) and Z ? (Z ? Z) (right). Observe that the peak is 
quite stable in the left-hand plot indicating the mean geodesic length is linear, while 
the right-hand plot the peak shows clear a left drift indicating that the geodesies 
grow sublinearly. 
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Figure 6: Plot of the mean geodesic length divided by n v \ for v = 0.98, 0.99 and 1. 
This data strongly suggests that Thompson's group has a non-trivial rate of escape. 
Note that the statistical error was smaller than the symbols used. 



Our estimates were quite sensitive to changes in 5: 



6 





1/4 


1/3 


1/2 


2 /3 


3/4 


A 


0.281 


0.279 


0.279 


0.276 


0.272 


0.267 


b 


176 


17 


8.0 


1.8 


0.47 


0.25 



Hence we conclude that the rate of escape is approximately 0.27 with an error of 
±0.01. 

We would like to conclude that this positive rate of escape implies that Thomp- 
son's group is non-amenable, however there are examples of amenable groups with 
non-trivial rate of escape. The group Z 3 I Z 2 is amenable but has positive rate of 
escape [20] . Unfortunately, computing geodesies in this group is equivalent to solving 
the traveling salesman problem on Z 3 [17] and so beyond these techniques. 



4 Conclusions 

We have computed exact lower bounds on the cogrowth of several groups including 
Thompson's group F. In particular, the cogrowth of Thompson's group must be 
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greater than 2.17329. By extrapolating the sequences of lower bounds we see that 
the bounds for the amenable groups clearly converge to 3, while those of the non- 
amenable groups converge to numbers strictly less than 3. Thompson's group appears 
to behave quite differently from the other groups we examined. Our extrapolations 
do not give clear results, though perhaps they point towards non-amenability. 

To further probe this group we used flat histogram methods to estimate the distri- 
bution of geodesic lengths in random words. The data suggests that geodesic lengths 
have an approximately Gaussian distribution about their mean length. Similar Gaus- 
sian distributions were observed for other groups, both amenable and non-amenable. 

The mean geodesic length of the amenable groups studied grow sublinearly, while 
those of Z I F 2 and Thompson's group are observed to grow linearly. Using simple 
sampling we estimate that the mean geodesic length of Thompson's group does indeed 
grow linearly and that the rate of escape is 0.27 ± 0.01. 
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